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Abstract 



We investigate the production of gluon pairs from a space-time dependent 
classical chromofield via vacuum polarization within the framework of the 
background field method of QCD. The investigation of the production of gluon 
pairs is important in the study of the evolution of the quark-gluon plasma in 
ultra-relativistic heavy- ion collisions at RHIC and LHC. 



I. INTRODUCTION 



Ultra-relativistic heavy-ion collisions at RHIC and LHC will provide the best opportunity 
to study the color deconfined state of matter, namely the quark-gluon plasma (QGP). The 
space-time evolution of the QGP can be split into different stages: 1. the pre-equilibrium, 
2. the equilibrium, and 3. the hadronization stage. One of the central problems in these 
experiments is to study how partons are formed and how their distribution function evolves in 
space-time to form an equilibrated quark-gluon plasma (if at all). High momentum partons 
(Pt > IGeV), i.e. minijets are calculated using pQCD. Soft Parton Production is treated 
differently. There exist various model approaches: 1) In the HIJING model soft parton 
production is treated via string formations. 2) In the color flux-tube model, an extension of 
the model named before, they are treated via the creation of a classical chromofield. When 
partons and a classical chromofield are simultanously present, a relativistic non-abelian 
transport equation has got to be solved. 



II. FIELD AND PARTICLE DYNAMICS 



The space-time evolution of the partons can be studied by solving relativistic non-abelian 
transport equations for quarks and gluons JI|,0. As the chromofield exchanges color with 
quarks and gluons, color is a dynamical quantity. The time evolution of the classical color 
charge follows Wong's equations M: 

^ = gf abc u,Q b A^. (1) 
There is also a non-abelian version of the Lorentz-force equation: 

^1 = g Q-F a ^u u (2) 



Taking the above equations into account, one finds the relativistic non-abelian transport 
equation 

faff* + gQ a F^ lu p v d^ + gf abc Q a Ayd c Q ]f(x,p, Q) = C + S. (3) 

Note that there are seperate transport equations for quarks, anti-quarks, and gluons. 
The single-particle distribution function f(x,p, Q) is defined in the 14-dimensional extended 
phase space of co-ordinate, momentum, and SU(3)-color. The first term on the LHS of 
Eq.(^|) corresponds to convective flow, the second to the non-abelian generalization of the 
effect of the Lorentz force, and the third term describes the precession of the color charge 
in the presence of a classical field. On the RHS there is the collision term C and the source 
term for particle production S. For any system containing field and particles one has the 
following conservation equation: 

d.TZt + W = (4) 

which is coupled with the above transport equation for the description of the QGP. 
The evolution of the plasma depends crucially on the source term S which contains all the 
information about how partons are produced from the classical chromofield. 



III. PARTON PRODUCTION FROM A SPACE-TIME DEPENDENT 

CHROMOFIELD 



The background field method of QCD is a suitable method to describe the production 
of partons from the QCD vacuum via vacuum polarization in the presence of a classical 
chromofield. Let us apply the background field method of QCD in order to describe the 
production of gg-pairs. 

The situation is similar to that of e + e~-pair production described by Schwinger [|J in 
QED. For a space-time dependent classical field A c i the amplitude for e + e~-pair production 
(see Fig. ([I])) from the vacuum is given by: 

M=<k li k 2 \S^\0> = -ieu(k 1 )^A^{K=k 1 + k 2 )v(k 2 ) (5) 
What, by the general formula: 



(2vr) 3 2A;0 (2tt) 3 2^ 
leads to the pair-production probability 0: 
r (i) a f ,a„„ 4m?, 5 _ 2m 2 



spin 



W^ e - = -Jd*K(l- (1 + -jj±)[\K ■ A d (K)\' - K'\MK)n (7) 

where the <i 4 iT-integral is defined for K 2 > Am 2 e . Simillarly, carrying out the same 
procedure in the non-abelian theory one finds for the amplitude for gg-pair production: 

M = iguXh^^Al^k, + k 2 )vi(k 2 ) (8) 

and for the corresponding probability |J: 

=Ey/^(l-^)V + ^)[l^ K(K)\ 2 - KM a d{K)\\ (9) 
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Fig. 1 Vacuum polarization 
diagram for the production 
of fermions in lowest order 



A. Gluon-Pair Production from a Space-Time Dependent Chromofield 

As conventional QCD cannot describe the interaction between a classical chromofield 
and a quantum gluon, one has to fall back on the background field method of QCD. This 
problem did not arise in QED, as there is no direct interaction between the classical field 
and the photon. That method was first introduced by DeWitt and further developped 
by 't Hooft ||. In the background field method of QCD, one defines: 

A afl = AJ + A**, (10) 

where A c i will not be quantized. So the generating functional excluding quarks is: 

Z[J, A d ] = J [dA q ] det Mq exp(i[S[A q + A d \ - ^G ■ G + J ■ A q \), (11) 
with the classical action: 

S[A q + A cl ] = -- A Jd 4 x(F a ^) 2 , (12) 
where the field-tensor is defined as: 

pa^u = g^ A au + A au^ _ qv ^ + ^ + g f abc ^ + A %)( A *> + ( 13 ) 

The gauge fixing term G a is chosen following 't Hooft: 

G a = Q^ A a^ + g jabc A ^ A ^ (14) 

The matrix element of Mq is given by: 

^•»» -=f Ssf (15) 

which is the functional derivative of the gauge fixing term with respect to the infinitesimal 
change of the gauge parameter 9 of the gauge transformation 
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5A an = _jabc gb^jw + A c^ + _ g^a 
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(16) 



Writing det Mq as functional integral over the ghost field, one obtains for the generating 
functional: 



Z[J,A d ,C,C]= j[dA q ][d X ]W 



x exp{i[S[A q + A cl ] + S ghost - —G -G + J ■ A q + X *£ + Cx}), 

la 



(17) 



where £ and are source functions for the ghosts and the ghost-part of the action is 
given by: 



S ghost = -J dW a p 2 S ab + g d ,r b \A c » + Af) 
-gf abc A^d, + g 2 f ace r db A^(A d c ? + Af)\ Xb . 



(18) 



Feynman rules involving a classical chromofield, gluons and ghosts can now be con- 
structed from the above generating functional ||. The vertices involving the coupling of 
two gluons to the classical field are given by: 



1A. 



\abd 
fiup 



gf abd [-2 gfip K u + g up (h - k 2 ) p + 2g^K p 



for coupling to the classical field once and by 

(V 2 aTX = -^ 2 [f abx f xcd (9,xg, P - g„g„x + g^gx P ) 



+f f (9^gx P - 9 P x9u P - 9 PP 9vx) 

+f acx f xbd (g^9x P - g„g v x)\ 



for coupling to the classical field twice. 



(19) 



(20) 
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Fig. 2 Vacuum polarization diagrams for the production of 
gluons in lowest order. 



Note that the above vertices are different from the three and four gluon vertices used in 
conventional QCD and that from hereon the classical field is denoted only by A not A c [. The 
gluon production amplitude M =< k\k2\S^ |0 > is defined in a way so that contains 
all interaction terms of the Lagrangian density involving two Q-fields, i.e.: 
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c(i) _ c(l) _i_ o(l) 

= iy d*x(-ii^[A]^/ ate g 6 »*g o ' 

™ - d v Q$gf*°(A*Q a ' + Q btJ, A cu ) 
-\g 2 r hc r h ' d {A h o c v + gM^)(A 6 '^Q cV + Q 6 'M cV )) 



1 d A x{-d x Q aX gf ahc A\Q CK 



--g 2 f abc f abc A b x Q^A b K Q CK ). (21) 

The total amplitude M = Ml4 + consists of a contribution by the three-vertex (see 
Fig-(@)(a)): 

M 1A = ^fjd*KS^(K-k 1 -k 2 ) 

A a »(K)e bl/ (ki)e d »(k 2 )(V 1A )f„ d p (22) 
and one by the four- vertex (see Fig.(|2|)(b)): 

M 2A = ^J (fk^hS^ih + k 2 -k 3 - k 4 ) 

A^ik^^k^ik^ih)^)^. (23) 

The above amplitudes include all the weight factors needed in order to retrieve the 
corresponding Lagrangian density. Now, we again calculate the pair production probability: 

-^JjSkwm^- (24) 

To obtain the correct physical gluon polarizations in the final state we use: 

E ^iK^O = E t u {k 2 y»\k 2 ) = -<r' (25) 

spin spin 

for the spin-sum and afterwards deduct the corresponding ghost contributions. 
The probability for the gluon part becomes: 

1 D r 

W 9 = —a s / d A K 
8 J 

((A a (K) ■ A* a {K))K 2 - (A a (K) ■ K){A* a {K) ■ K)) 

f aa ' c '[(A a (K) ■ A* a '(-k 3 ))(A* c '(K - k 3 ) ■ K)) 

+ ^|! J d A hd%d A K 
((A a (k 3 ) ■ A C (K - h))(A™'(k' 3 ) ■ A* C \K - k' 3 )) 

^jabx jxcd _|_ jadx jxcb^ja'bx' jx'c'd _|_ j-a'dx' jx'c'b^ 

(A a (k 3 ) ■ A* a \k' 3 ))(A c (K - k 3 ) ■ A* C '(K - k' 3 ))). (26) 
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Fig. 3 Vacuum polarization diagram for the production of 
ghosts in lowest order 



Now, we calculate the ghost part. The vertices involving two ghosts and one classical 
field and two ghosts and two classical fields respectively are given by: 

(Yf/Y™ = +9f abd (h - k 2 )» (27) 

and: 

(\rFP\abcd „'„2„ / rabx fXcd i fadx rxcb\ {OQ\ 

K V 2A = l 9 9fj,x{J J +J J )■ (/°) 
The corresponding amplitude for the ghosts reads: 

(M FP ) bd = {M[/) bd + (M£f) M (29) 

with (see Fig. ©(a)): 

(M[/) bd = J d 4 K(2n) 4 5' 4 \k 1 + k 2 -K) A a »(K)(V FP )f d (30) 

and (see Fig. ©(b)): 

(M FP ) hd = i J d 4 k 3 d 4 k 4 5 {4 \k 1 + k 2 -k 3 - fc 4 ) 

A a ^{k 3 )A c \k 4 ){W FP ) a ^ d . (31) 
The probability in this case is simply: 

tFP f d 3 ki d?k 2 / T . /r FP\bd/ 1 \ / rFP\*bd 



WFF - 1 jm4 (M ) iM r (32) 



which becomes: 

rFP _ a s 



W 1 



J d 4 K5 {4 \K-k 1 - k 2 ) 



((A a (K) ■ A a (K))K 2 - (A a (K) ■ K)(A a (K) ■ K)) 

asg2 f d 4 Kd%d% 



32 

(A a (k 3 ) ■ A\K - k 3 ))(A a '(k' 3 ) ■ A C '(K - k' 3 )) x 

^pbx pcd _|_ jadx jxcb^^jta'bx! jx'dd | ja'dx' jx'c'b^ (33) 
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The real gluon-pair production probability is given by W gg = W 9 — W FP . 

Instead of the probabilities for pair production, one can also consider the corresponding 
source terms which then ultimatively enter the transport equation. The source terms are 
equal to the probability per unit of time and per unit volume of the phase space. Some 
calculations yield |fT0|1 : 



dW$ g 2 m 
d 4 xd?k (27r) 5 o; 



"" ~ '' "' A°(x) e lk ' x J r/'r, A a u {x 2 ) e~ ih ' 



{i[W{x - x 2 y + {x- x 2 fk v + k ■ (x - x 2 )g^} 
,K (m\J—(x - x 2 ) 2 )m^J -{x - x 2 ) 2 + 2K\{m^j -{x - x 2 ) 2 ) 



KAmJ — (x — x 2 ) 2 ) 

•™y y, * - )■ (34) 

sj-{x - x 2 ) 2 



for the quarks and: 



dW gg 1 



l\2 



d 4 xd 3 k (27r) 5 fc° J (x - x') 

x{~g 2 A a »(x)A a »'(x')[3k^ - 8g^k vA ' 



4 [x — x 



A2 



.{x-x')^ .(x-x')^ 6g^ (x - £%(x - x')^ 



x — x') 2 (x — x') 2 (x — x') 2 (x — x') A 

-3ig 3 A a ^{x')A cX {x')A a ^'(x)f a ' ac K x g^ 

- g 4 A a11 (x)A cX (x)A a '»' (x')A c ' x ' {x') 

r_ „ / rabx rxcd , radx rxcb\( ra'bx' rx'c'd , ra'dx' rx'c'b\ 

[^aVa'U / + J J ){J J +J J > 

+^g^gxx>F cx f a ' c ' x }}. (35) 

for the gluons. It can be checked that the above results are gauge invariant with respect 
to type-(I)-gauge transformations [ffO . 



IV. DISCUSSION 

The above results are still to complicated in order to directly get an idea about their 
content, so we look at them for a special, purely time dependent model field. 

A a \t) = A m e-W /t0 , t > 0, a = 1, ...,8, (36) 

and all other components are equal to zero. Many other forms could have been taken. 
We have chosen this option just to get a feeling for how the source term in the phase-space 
behaves. The actual form of the decay of the classical field can only be determined from a 
self consistent solution of the relativistic non-abelian transport equations. The above choice 
yields: 
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Fig. 4 Source term S [MeV] for quarks (dashed) and gluons (solid) production respectively versus 
time t [fm/c] and time- integrated source-tem T for quarks and gluons versus kx [MeV] for the 
above choice of the model field. 



= ^(^V^-N*-* (37) 



d A xd 3 k (2tt) 2 v m ' 1 + AuHl uj 2 ' 

with m 2 ^ = m 2 + k\ where kr is the transverse momentum and: 

dW 99 2Aa s 2 2ik o t lt{/to t fe| 
d*xd 3 k (2vr) 21 m) l + 4(k°)H 2 { (P) 2J 

36a| 4 2ikH 2 |t|/to *o 1 /oon 

+ 2vr (Anj 6 6 l + (Jfc0)^(Jfc0)a- ^ 

We choose the following parameters: as = 0.15, ^4 in = 1.5GeV, k T = 1.5GeV, y = 0, 
and t — 0.5/m. Additionally, the quarks are considered to be massless. On the LHS of 
Fig.([|), the oscillatory behavior of the source terms S seems to indicate that there exist 
periods of particle creation and particle annihilation which follow each other periodically. 
This oscillatory behavior of the source term will play a crucial role once it is included in a self 
consistent transport calculation. It can also be seen in the Figure that there are considerably 
more gluons produced than quarks. On the RHS of Fig.(^), the time-integrated source terms 
T can be regarded as a measure for the net-production of particles in an infinitesimal volume 
around any given point in the phase-space. It does not show the oscillatory behavior which 
gives a totally different picture for different times. In future, we will include these source 
terms in the transport equation in order to study the production and equilibration of the 
QGP at RHIC and LHC. 
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